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It is perhaps appropriate that, in a year marking the 90th anniversary of Meghnad 
Saha seminal paper (1920), new developments should call fresh attention to the 
problem of ionization equilibrium in gases. Ionization equilibrium is considered in 
the simplest "physical" model for an electronic subsystem of matter in a rarefied 
state, consisting of one localized electronic state in each nucleus and delocalized 
electronic states considered as free ones. It is shown that, despite the qualitative 
agreement, there is a significant quantitative difference from the results of applying 
the Saha formula to the degree of ionization. This is caused by the fact that the 
Saha formula corresponds to the "chemical" model of matter. 

PACS number(s): 52.25.Jm, 52.25.Ya, 05.30.Fk, 34.80.Dp 

The Saha equation [1] derived in 1920 plays a fundamental role in all fields of plasma 
physics [2], and defines the ratio of electrons, ions, and atoms at given thermodynamic 
parameters in an ideal (or weakly nonideal) mixtures of these particles. The Saha equation 
also plays an important role in cosmology, defining the recombination transition in models 
of the early Universe [3,4]. 

The Saha equation was derived based on the matter model in which bound states of 
particles are considered as new composite particles with a priory specified structure and 
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properties. This means, that real matter which should be considered in most cases as a 
quasi-neutral non-relativistic many-body system consisting of interacting nuclei and elec- 
trons ("physical" model of matter [5]) can be considered under certain conditions as a 
system of atoms or molecules. Furthermore, atoms or molecules in such a model, being in 
essence bound states of a certain finite number of electrons and nuclei, represent a form 
of quasiparticles [5]. A statistical description under such consideration requires a strong 
assumption which is a basis of the so-called "chemical" model of matter, i.e., fundamental 
(electrons and nuclei) and composite particles should be considered equivalently (see, e.g., 
[5]). Attempts to justify this assumption are not terminated to date (see, e.g., [6] and ref- 
erences therein). In our opinion, this is caused by two important factors. One situation is 
associated with that the " chemical" model has reasonable physical grounds when considering 
matter in the rarefied state. However, in describing matter at high densities, the concept on 
atoms, molecules, and other composite particles loses meaning. This means that the statis- 
tical description of the " chemical" model should imply the formalism of the " appearance" 
or "disappearance" of composite particles (atoms and molecules). In this case, the case in 
point is their principal presence or absence, rather than their small number which is assumed 
by the Saha formula. Another situation is associated with the consideration of the iden- 
tity (indistinguishability) of fundamental particles, first of all, electrons (see, e.g., [7]), whose 
consequence is the Pauli principle for electrons [8] . The point is that the number of chemical 
potentials in the statistical description of matter in the "physical" model (which is primary) 
is defined by the number of types of distinguishable particles (see, e.g., [9]). In particular, 
for the pure matter consisting of electrons and nuclei of the same type, the theory implies 
two chemical potentials, i.e., the chemical potential of electrons and the chemical potential 
of nuclei. In this case, three chemical potentials are put into consideration in the "chemical" 
model, i.e., the chemical potential of "free" (delocalized) electrons, the chemical potential of 
" ions" (or free nuclei) , and the chemical potential of " atoms" consisting of bound electrons 
and nuclei. Thus, if the problem of the introduction of the chemical potential of "atoms" 
(see above) is not considered, the possibility of replacing the initial chemical potential for all 
electrons by the chemical potential of only "free" electrons remains unclear. In our opinion, 
the only possible justification of such a step from the viewpoint of indistinguishability of 
electrons is the recognition of the fact that, although electrons themselves are indistinguish- 
able, electronic states are distinguishable (see, e.g., [10]). In particular, this is vahd for the 
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distinguishability of localized ("atomic", "molecular", etc.) and delocalized ("free") elec- 
tronic states. This circumstance formally offers the possibility of "different" mathematical 
descriptions of localized and delocalized electronic states using various chemical potentials. 
From this point of view, the actual reason of the development of the "chemical" model of 
matter is caused by difficult construction of the general theory of systems with Coulomb 
interaction. The point is that it is necessary to uniformly describe both localized electronic 
states which are characterized by the strong electron-nucleus interaction and delocalized 
electronic states which are quite adequately described within the perturbation theory with 
respect to interparticle interaction (see, e.g., [9]). In this situation, the problem of the quan- 
titative relation between the results of applying the "chemical" and "physical" models to 
describe the ionization equilibrium in the rarefied gas state, where the "chemical" model 
has reasonable grounds, becomes central. In the case of the "physical" model, we shall 
proceed from the assumption on the classical description of the nucleus subsystem, which 
corresponds to the Saha formula derivation. The electron identity principle requires that 
all electrons of the system should be described by the uniform chemical potential He- When 
considering localized electronic states in the low-density limit, we restrict the analysis to 
the one-center approximation which is a necessary condition of the existence of " atoms" as 
quasiparticles. For convergence of statistical sums, the electron-electron interaction can be 
considered within the self-consistent Hartree-Fock approximation [11]. 

Denoting electron energy levels in locahzed and delocahzed states by En and €{q), re- 
spectively, and assuming that localized states are take place at each nucleus, for the average 
number {N^Y^^'' of electrons in the system, we obtain [12] 



Here Nc is the full number of nuclei in the system. The function fe{E) = [1 + exp(£^ — 
lie)/T]~^ is the Fermi energy distribution. The first and second terms on the right-hand side 
of (1) yield the average number of electrons in localized {N^'^^°^'> and delocalized (iVg) 
states. In this case, Eq. (1) together with the quasineutrality condition (Ng)^'^^^ — ZcNc 
makes it possible, for a given spectrum of electronic states, to determine the chemical po- 
tential of electrons as a function of the nuclei density Nc/V, temperature T, and nucleus 
charge Zc, hence, to determine the degree of ionization a equal to ratio of the number of 




(1) 
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delocalized electrons to the number of nuclei in the system under study, 



a — 



(2) 



Bearing in mind the comparison of this approach with the Saha formula for the ionization 
equilibrium in its simplest form, we make a number of simplifying (but generally speaking, 
unnecessary for the general approach described above) approximations. Let us assume that 
delocalized states in the low-density limit can be approximately described by plane waves 
with e(g) = h?q^ /2me (strictly speaking, delocalized states orthogonal to localized electronic 
states should be searched, e.g., in the form of the so-called COPW states [13]). Furthermore, 
when considering delocalized electronic states, we suppose that /Xg < 0, | [Xe/T \» 1, 
replacing the Fermi distribution by the Boltzmann distribution, whereas localized states 
are certainly described by the Fermi distribution. In this case, the number of delocalized 
electronic states in the volume V is 2V exp(/Xe/T'). Equation (1) for determining the chemical 
potential is still transcendental, and the problem of its numerical solution is related to the 
necessity of determining the energies of localized excited (n > 1) electronic states in the 
Hartrcc-Fock approximation [12]. 

We will now use the known result based on the solution of the Saha equation. When 
determining the degree of ionization for hydrogen and of some other elements, this result 
makes it possible to restrict the analysis to the consideration of only the electron ground 
level in the atom and states of the continuous spectrum in wide temperature and density 
ranges, disregarding excited localized states [14]. In this case, putting the excited states into 
the calculation slightly affects the accuracy of the determination of the degree of ionization 



Let us perform the further consideration without loss of generality; for the case of hy- 
drogen, Zc — l- Under given conditions, the ground energy level for the localized electronic 
state £'o = — / in the Hartree-Fock approximation coincides with that for the case of com- 



to the Bohr atom) [12]. In this case, it follows from (1) that the chemical potential in the 
approximation under consideration is explicitly defined. In the case of hydrogen, = 1, we 



[15]. 



plctc disregard of the electron-electron interaction, when / = mgC' 



/2/?.^ (which corresponds 
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obtain the equation for //g 



ii, -T\^\\ -[exp(-(7/T) + ^] + y [exp(-//T) + •^\^ + 2neA3exp(-7/r) \ ; 



Based on (3), the degree of ionization is calculated as 



^ -[1 + ^ exp(//T)] + ^[1 + !-^i£^]2 + 2neA3 exp(//r) 
exp(-7/T) L_ (4) 



In deriving (2) and (3), we considered the spin degeneracy factor for energy levels (which 
was ignored in [12]). 

As is known, the degree of ionization calculated in the same approximation by the Saha 
formula [1,14], under the assumption of ideality of atomic, electronic, and ionic components 
of ionized gas, is given by 

a(^) = 2 exp(-//T) + + ^"^^^ (5) 

As it is easy to see, both formulas for the degree of ionization are identical in the low-density 
limit under the condition that value 7 = ncA3exp(//T) ^ 1. In this case and at a fixed 
temperature, the degree of ionization is close to full ionization with the accuracy to first 
order on the value 7, 



lim a ~ ~ 1. (6) 

The same result corresponds to the transition to high temperatures at a fixed density of 
nuclei. 

When passing to low temperatures, at a fixed density, the degree of ionization is close to 
zero, 

lim a ~ q;(^) = 0. (7) 

T->-0,n 
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Thus, it is clear that the point n = 0, T = is a singular point at the degree of ionization a, 
similarly to that the point ^ = 0, a; = is a singular point for the permittivity e{q, uj) which 
depends on the wave vector q and frequency u (see, e.g., [16,17]). This means that matter 
will be in an atomic state at a fixed low density n in the limit T — )■ 0. In turn, at a fixed 
temperature T in the limit n ^ 0, matter will be in a completely ionized state. Although 
the points n — and T — are practically inaccessible, the non-permutability of limits (see 
also [15]), written for the degree of ionization as 

lim a{nc,T)j^ lim a{nc,T). (8) 

T->0,nc->-0 nc-*-0,T->0 

has important physically observable consequences. 

Figures 1 and 2 show the degrees of ionization and chemical potentials in the physical 
model (solid lines) and according to the Saha formula (dashed lines) for identical approxi- 
mations (the consideration of only one, i.e., the lowest, localized state and the consideration 
of delocalized electrons as free particles) . 

We note that since ions are originally absent in the physical model, their average charge, 
can be additionally defined either as the number of free electrons per one nucleus not occu- 
pied by electrons. 



or as the number of free electrons relative to the total number of nuclei. 



^^ - ^V^' (10) 

For hydrogen Zi — 1, while z* — a and varies from to 1, being a new statistical variable in 
the electron-nucleus system. 

As follows from the above calculations, despite the qualitative agreement of the results 
of applying the "chemical" and "physical" models to the degree of ionization in rarefied 
matter, there is a rather serious quantitative difference between them even in the considered 
simplest case. The advantage of the physical model consists in enumeration of identity of 
all electrons in the system. In fact, in the physical model, electrons are distributed on all 
existing electron states, according to the Pauli principe, in contrast with the chemical model. 
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where electrons at first are separated on atomic and delocafized states, which are described 
by a different way. It is easy to estabhsh a formal similarity between application of the 
physical model to plasma ionization in the present paper and consideration of electrons in 
semiconductors with donor and acceptor impurities (sec, e.g., [18]). 

The approach firstly used in the present paper can lead to very significant quantitative 
differences in the calculation of thermodynamic, kinetic, and electromagnetic properties, 
since the degree of ionization (or the number of delocalized ("free") electrons) is of funda- 
mental importance. Thus, the consistent development of the "physical" (Coulomb) model 
is of paramount problem for studying real matter even in the rarefied state. 
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Figure 1: Ionization a for hydrogen as function of dimensionless temperature x = T /I for the value of the 
parameter tIcAq — 0.01. The variable Aq = (27r?i^/m/)^/^. Solid line - physical model, dashed line - the Saha 
formulae. 
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Figure 2: Dimensionless chemical potential /ie// for hydrogen as function of dimensionless temperature 
X = T/I for the value of the parameter ricAg = 0.001. Solid line - for all electrons in the physical model, 
dashed line - for free electrons, accrording to the Saha formulae. 



